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Abstract: We carry out an exact analysis of the average frequency v^x in ^ ne direction Xj of positive- 
slope crossing of a given level a such that, /i(x, t) — h = a, of growing surfaces in spatial dimension d. 
Here, /i(x, t) is the surface height at time t, and h is its mean value. We analyze the problem when 
the surface growth dynamics is governed by the Kardar-Parisi-Zhang (KPZ) equation without surface 
tension, in the time regime prior to appearance of cusp singularities (sharp valleys), as well as in the 
random deposition (RD) model. The total number N + of such level-crossings with positive slope in all 
the directions is then shown to scale with time as t d l 2 for both the KPZ equation and the RD model. 
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1. Introduction 



Due to their practical applications and fundamental interest, a great amount of effort has been 
devoted to understanding the mechanism(s) of growth of thin films, and the kinetic roughening of 
their surface during the growth. Analytical as well as numerical analyses of many models of surface 
growth suggest, quite generally, that certain properties of the surface exhibit dynamic scaling. 1-6 In 
order to derive quantitative information about the films' surface morphology, one considers a sample 
of size L and defines the mean height h of the growing film and its surface width w by, 1 

1 r L / 2 

h(L,t) = - dx h(x,t) , (1) 

L J-L/2 

where h(x, t) is the surface height, and 

w(L,t) = [((h-h) 2 )] 1/2 , (2) 

with (•) indicating an average over different realizations of the surface. Starting from a flat surface 
as the initial condition, it was proposed by Family and Vicsek 7 that rescaling the space and time 
variables by, respectively, b and b z rescales the surface width w by b x , w(bL, b z t) = b x w(L, t), implying 
that 

w(L,t)=L*f(t/L z ), (3) 

where f(x) is a scaling function, and z is the dynamic exponent. If for large t and fixed L (t/L z — > oo) 
the width w saturates, then one must have, f(x) — ► constant as x — ► oo. However, for fixed and large L 
and 1 <C t <C L z , one expects the correlations in the height fluctuations to grow within only a distance 
t l l z and, thus, w must be independent of L. This implies that for x <C 1, f(x) ~ x 13 with, [5 = %/z. 
Therefore, the dynamic scaling of Family and Vicsek 7 postulates that, w(L,t) ~ t@ for 1 ^ t <S L z , 
and w ~ L x for t 3> L z . The roughness exponent \ an d the dynamic exponent z characterize the 
self-affine geometry of the surface and its dynamics, respectively. These scaling relations have been 
tested for a variety of growing surfaces by extensive numerical simulations and analytical calculations. 

In this paper we introduce the concept of level crossing in the context of surface growth processes. 
In the level-crossing analysis one is interested in determining the average frequency, v^xi in the 
Xi— direction of observing a given level a for the function, h — h = a, in the growing film. We show 
that v+ x . is written in terms of the joint probability distribution function (PDF) of h — h and the 
gradient of h(x, t). Therefore, v+ x . carries the same information about the film's surface which is 
contained in the joint PDF of the fluctuations in the height /i(x, t) and its gradient. Our goal in this 
paper is to study and analyze the quantity v£ x . at time t for a growing surface in a sample of size L. 
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In addition, we introduce an integrated quantity, iV+, defined by, iV+ = J^da v+ x ., which mea- 
sures the total number of positive-slope crossings of the surface in the Xj— direction, and is expected to 
become size-dependent in the stationary state. We determine exactly the time- and height-dependence 
of v&xi an d N£. for the Kardar-Parisi-Zhang (KPZ) equation in the limit of strong coupling, over time 
scales prior to the emergence of cusp singularities (sharp valleys). We also derive exact expressions 
for the same quantities for the random deposition model. 

The statistics of level-crossings that we study in the present paper are closely related to the concepts 
of first passage and persistence of fluctuating interfaces in both space and time, which have recently 
been studied in many papers. 8,9 In particular, the concept of spatial persistence was introduced by 
Majumdar and Bray 10 (MB), who showed that the probability Pq(1) that the height of a fluctuating 
(d + 1)— dimensional surface in its steady state stays above its initial value up to a distance I, along 
any linear cut in the d— dimensional space, decays as, Pq(1) ~ l~ 8 ■ Here, 9 is a spatial persistence 
exponent which takes on distinct values, depending on how the point from which the distance I is 
measured is specified. Their analysis was carried out for a fluctuating interface at steady state. The 
exponents were shown to map onto the corresponding temporal persistence exponents for a generalized 
(d + 1) — dimensional random walk. 

There are a few important differences between our work and that of MB. First, all the analysis 
and calculations that we present in the present paper are valid at times t < t*, where t* is the time 
scale over which sharp valleys (cusp singularities) are developed in the growing surface. Therefore, the 
time regime that we investigate in the present paper is completely different from what was considered 
by MB, as we study the probability of crossing a certain height with a positive slope at times t < t* , 
for which we study the short-time behavior of the KPZ equation. Secondly, we study the KPZ model 
in the limit of zero surface tension - a nonlinear equation - and derive our exact results for any spatial 
dimension d, whereas MB derived their results for a linear model (the Gaussian model) , except when 
they considered the KPZ equation in (l+l)-dimensions. Thirdly, we study the problem assuming a 
flat initial surface, which is distinct from "finite initial starting point" of MB (see below). 

The plan of this paper is as follows. In Section 2 we discuss the connection between u^ x . and the 
underlying PDF of growing surfaces. In Section 3 we derive an integral representation for v~^ x . for 
the KPZ equation in (d + 1)— dimensions in the strong coupling limit, before the emergence of the 
cusp singularities. Section 4 presents our exact results for the random deposition model. The paper 
is summarized in Section 5, while four Appendices provide the details of the derivation of our results. 
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2. Level-Crossing in Growing Surfaces 



Consider a sample of an ensemble of functions which make up the homogeneous random process 
/i(x, t), representing the height of a growing surface. Let n+ denote the number of positive-slope 
crossings such that, /i(x, t) — h = a. In a time t and for a typical growing surface of linear size L (see 
Figure 1), let the mean value of n+ for all the samples be N+(L), 

N+(L) = E[n+(L)], (4) 

where E denotes the expectation or mean value of the quantity. If we consider iV+ in a second segment 
of size L immediately following the first, then, since the process is homogeneous, we obtain the same 
result as given by Eq. (4). Thus, for the two intervals together we obtain, 

N+(2L) = 2N+(L) , (5) 

from which it follows that, for a homogeneous process, the average number of crossings is proportional 
to the space interval L. Hence, 

N+(L) oc L , (6) 

or 

N+(L)=v+L. (7) 

We now consider how f + is deduced from the underlying probability distribution for /i(x, t) — h. 
Consider a small length dl of a typical sample function. Since we assume that the process /i(x, t) — h 
is a smooth function of x with no sudden increase or decrease, and that dl is small enough, then the 
sample can only cross the level /i(x, t) — h = a with a positive slope, if h(x, t) — h < a at the beginning 
of the interval dl. Furthermore, there is a minimum slope at position x if the level h(x.,i) — h = ais 
to be crossed in the interval dl, depending on the value of h(x,i) — h at x. Therefore, there will be a 
positive crossing of /i(x, t) — h = a in the next space interval dl, if at position x, 

h(x,t)-h<a, and > <* = ^ j) Z S , (8) 

If the above conditions are satisfied, then there will be a high probability of crossing the level in the 
interval dl. 11 ' 12 

In order to determine whether conditions (8) are satisfied at an arbitrary location x, we must 
determine how values of y = h(x,t) — h and y' = dy/dl are distributed by considering their joint 
probability density P(y, y'). Suppose that the level y = a and the interval dl are specified. Then, we 
are interested only in those values of y such that y < a and of y' = dy/dl > (a — y)/dl, which represent 



the region between the lines y = a and y' = (a — y)/dl in the plane (y,y'). Hence, the probability of 
a positive-slope crossing of y = a in the interval dl is given by, 

roo ra 

/ dy' dyP(y,y'). (9) 

JO Ja-y'dl 

As dl — ► 0, one has, 

P(y,y') = P(y = a,y') . (10) 
Since for large values of y and y' the PDF approaches zero fast enough, Eq. (9) may be written as 

roo pol 

\ dy' dyP(y = a,y'), (11) 

JO Ja-y'dl 

in which the integrand is no longer a function of y, so that the first integral in Eq. (11) is simply, 
Ja-y'dl dy P{y = a >y') = P(y = a -,y')y'dl- Then, the probability of a positive-slope crossing of y = a 
in the interval dl is equal to 

roo 

dl P(a,y')y f dy' . (12) 
Jo 

Since according to Eq. (7) the average number of positive-slope crossings over a length scale L is 
f^L, then, the average number of crossings in the interval dl is v^dl. It then follows that the average 
number of positive-slope crossings of the level y = a in the interval dl is equal to the probability of 
positive-slope crossing of the level y = a in dl, which is true only if dl is small and the process y{x) is 
smooth enough that there cannot be more than one crossing of y = a in the interval dl. In that case, 
we have v+dl = dl / °° P(a, y')y'dy' and, therefore, 

P{a,y')y'dy' . (13) 



In the following sections we will derive exact expressions for via the joint PDF of /i(x, t) — h 
and the height gradient, for both the KPZ and random deposition models. To derive the joint PDF 
we use the master equation method. 13-16 This approach enables us to determine m terms of a 
generating function. For example, the generating function for a (2 + 1)— dimensional surface is given 
by (see below), 

Z(X, /j,, x, t) = (exp{— iA[/i(x, t) — h] — ^^(x, i)}) , (14) 

where, it(x, t) = — V/i. 

3. Analysis of Level- Crossing for the KPZ Equation 

In the KPZ model in (d + 1)— dimensions, the surface height /i(x, t) at position x on top of the 
substrate, in the limit of the zero surface tension, satisfies the following stochastic equations 

f t =h t =\*i: «?+/, as) 
i=i 



-^1 = u i:t = a^2 u iPji + fxi , (16) 
i=i 

where u% = dh/dxi = h Xi , and pij = dh x Jdxj. Here, / is a zero-mean random force with a Gaussian 
correlation in space and white noise in time: 

</(x, t)/(x', 0) = 2£>„£>(x - x!)6(t - t') , (17) 

where D(x — x'), an even function of its argument, is the spatial correlation function which takes on 
the following form, 

D <*-*'> = ^ I .... exp[-£fei-^g 

with a Xi being the standard deviations in the Xi— direction. The parameters a and Dq describe, 
respectively, lateral growth and the noise strength. Typically, in order to account for short-range 
correlations, the correlation function is taken to be, Z)(x — x') = <5(x — x'), but we regularize the 
delta-function correlation by a Gaussian function. When the standard deviations a = a Xi are much 
smaller than the system's size L, we would expect the model to exhibit short-range correlations in the 
/ term. Therefore, we emphasize that our analysis is for the case when, a <^ L. 

As analyzed in Appendices A, B, and C, we define the generating function, Z(X, fii, Xi,t) = 
(0(A, Hi, Xi, t)), for the fields h = h(x,t) — h and Ui = h Xi , where 



(18) 



= exp |— iX [/t(x, t) — h(t)] — • 



(19) 



A and \ii are the sources of h and u^, respectively, and, i,j = 1, ■ ■ ■ , d. Assuming statistical homogene- 
ity, i.e., assuming that, dZ/dn. = 0, it follows from Eqs. (15) and (16) that Z satisfies the following 
equation (see Appendix C for details), 

BZ 1 

Zt = ~dt = lXl{t)Z ~ 2 lX ^ Z ^ 1 " x2k (°) Z zZrik"(0)Z . (20) 

where, fc(x-x') = 2D D(x-x'), 7 (t) = h t , k(0) = 2D /(ira d ), and ^.(0) = -4D /(ira 2+d ), where 
we use, a = a Xi = a Xj , for simplicity. 

In trying to develop a statistical theory of level-crossings in rough surfaces, it becomes clear that 
the interdependence of the statistics for height difference /i(x, t) — h and height gradient must be taken 
into account. The very existence of a nonlinear term in the KPZ equation leads to development of 
the cusp singularities (sharp valleys) in a finite time and in the strong coupling limit, hence forcing 
one to distinguish between different time regimes. It was shown recently that, starting from a flat 



6 



surface, the KPZ equation will develop sharp-valley singularities after a time scale t* where, 13 t* ~ 
Do ~ 1 / 3 a~ 2 /'V( d+4 )/ 3 . This implies that for times t < t* the relaxation contributions tend to vanish 
in the strong coupling limit. In this regime one can derive closed- form solution for the generating 
function. Starting from a flat surface, i.e., from h(pc, 0) = and it(x, 0) = 0, one has the following 
solution (see Appendix C) 



Z(X,m,---,H d ,t) = F^X, Ln,t) ■ ■ ■ F d (\, fi d ,t) exp[-\ 2 k(0)t} , 



(21) 



with 



i(X,Hj,t) = | 



1 - tanh 2 



yj2ik xx (0)a\t 



1. 2 2ik xx {0) 

x exp < iUj\ — ; tanh 

2 J V aX 



\J 2ik xx (0)aXt 



1 



-1/4 



icek xx (0)Xt 2 



(22) 



To derive a closed expression for the PDF P(h,Ui,t), one needs to determine such moments as 
(h n u r i n u l jPij). As shown in Appendix B, such moments are identically zero. 14 Using this result, it can 
then be shown that P(h,Ui,t) takes on the following expression (see Appendix C), 

P(h,u Xi ,t) = ^y+i J dXdm ■ ■ ■ d/j, d exp ^iXh + i ^ HiU^j Z(X, m, ■ ■ ■ n d , t) . 



(23) 



According to the Eq. (13), the frequency of crossing a definite height h(x,t) — h = a with a positive 
slope in, for example, the x\— direction, is given by, 



/*oo /*oo 

v txi = / du *i u *i / du x 2 ■ ■ ■ du Xd P(a, u Xl , • ■ ■ , u Xl , t) . 

JO J — oo 



(24) 



Using the Eqs. (21)-(23), v£ x is then given by, 



ax i 



= — r 



1 

oo 



oo /-oo 



iXa 



OO J — oo 



Z(A, /ii , -»■ 0, ■ ■ ■ , // d -»■ 0, t)dXdm , 



dAexp iAa — X 2 k(0)t — iak xx (0)Xt 2 (, 



(25) 



where 




'1 - tanh 2 



tJ2ik xx (0)aX 



(26) 



To compute the integral in Eq. (25), we used the numerical integration software developed by 
Piessens et al. 17 which uses a globally-adaptive scheme based on Gauss-Kronrod quadrature rules. 
In Fig. 2 we plot i/+ Xl for times, t/t* = 0.05 ,0.1, and 0.2, which are before the emergence of the 
cusp singularities. To derive an expression for A+ , the total number of level-crossings with positive 



slopes in the x\— directions, let us express it in terms of the generating function Z. It can be shown 

straightforwardly that iV+ is written in terms of the generating function Z as 

1 f°° 1 C 
N+=— 2 Z ( A ^°'/ u i^2^0,---, W ^0,t)d /[ /i = lirn^ . (27) 

Using the Eq. (26) one finds that, iV+ ~ i 1 / 2 . In Fig. 3 (obtained by direct numerical integration of 
-^xi = X^oo v ax\ da) we plot iV+ vs t, which also indicates that, iV+ ~ i 1 / 2 , in agreement with the 
analytical prediction. 

4. Analysis of Level- Crossing for the Random Deposition Model 

In the random deposition model, the height of each column performs an independent random 
walk. This model leads to unrealistically rough surfaces, with its width growing with the exponent 
(3 = x/ z = 1/2 without ever saturating. In the continuum limit the random deposition model is 
described by, 

d d 

— h(x,t) = /(x,t) , —Ui(x,t) = f Xi , i = l, •••,<*, (28) 

where, Uj(x, t) = dh{x,t)/dxi, and /(x, t) is a zero-mean random force described by Eqs. (17) and 
(18). Assuming homogeneity, we define the generating function by 



Z(X,fj,i,t) = ( exp 



-i\h(x.,i) -i^mui 



(29) 



i=i 

As before, A and ^ are the sources of h and Ui, respectively, and i, j = 1, • • • , d. 

It is then straightforward to derive the following equation for the evolution of Z(X, Hi,t): 
f> ^ 

-Z(X, t) = -X 2 D D(0)Z + ]T n 2 iD Q D x . x . (0)Z . (30) 

i=l 

The joint PDF of h and ui is obtained by Fourier transform of the generating function: 

P(h,u,t) = — f dXdfii exp ( iXh + i \XiUi ] Z(X,fii,t) . (31) 
27rJ \ f^i J 

Carrying out the Fourier transformation, we obtain a Fokker-Planck (FP) equation, 

§P = D „ D (0)J^P-£lVW0)|jP. (32) 

2=1 1 

The solution of the above FP equation is written as, P(h,Ui,t) = pi(h, t)p2{u\, t) ■ ■ ■ Pd+i(ud, t)(for 
motivation see 20 ). Using the initial conditions that, pi(h, 0) = 5(h), and, Pi(v,i, 0) = 5(u{), and starting 
from a flat surface, it can then be shown that, 



where, D" = D XiXi = D XjXj . Thus, the frequency of crossing a definite height, fo(x, t) = a, in, for 
example, the x— direction, is given by, 



roo 

/ u x P(a,u x ,Uj)di 
Jo 

rv 2 ! n/ 2 

, j = 2,---,d. (34) 



D xx (0) 

exp 



2irV D{0) 



u 

1 



4D(0)t 



27TO- 



exp 



4D D(0)t 



Thus, the quantity v^ x in the random deposition model has a Gaussian form with respect to a. 
Moreover, it is easily seen that iV+ ~ i 1 / 2 , in agreement with the result for the KPZ equation. 

5. Summary 

In this paper we derived exact results for the statistics of level-crossing with positive slopes for 
growing surfaces that are governed by the KPZ equation in (d+1)— dimensions with a Gaussian forcing 
term which is S— correlated in time and contains short-range spatial correlations, in the limit of zero 
surface tension. The integral representation of the frequency of crossing, v^., in the Xj— direction is 
given for the KPZ equation in the strong coupling limit before the emergence of sharp valleys (cusp 
singularities). We also determined the quantity iV+ = f*™ da v+ xo which measures the total number 
of positive-slope crossing of the growing surface in the x«— direction, and showed that, iV+ ~ t 1 / 2 . 
Using statistical homogeneity, it is then clear that, v+ x . = v^ x p an d that, N + = iV+ • • • iV+ ~ t d l 2 . 
We also derived exact expressions for similar properties in the random deposition model. 

The ideas and techniques presented in this paper are quite general and may be used to determine 
v+ x . for processes that are governed by a general Langevin equation with given drift and diffusion 
coefficients. 
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ht = ±aJ2 u l + f, ( 37 ) 
=1 



Appendix A 

In this and the following Appendices we derive the governing equation for the generating function, 
Eq. (20). In particular, in the present Appendix we investigate a more general definition of the 
generating function in (d + l)-dimensions, and will further expand the study in Appendix B in order 
to derive an identity which is crucial for deriving Eq. (20) in Appendix C. 

We define a generating function by 

Z(X, m, r)ij,Xi, t) = (0(A, fj,i,rjij,Xi,t)} , (35) 

for the fields, h = h — h, Ui = h Xi , and pij = h XiXj . A, fii, and rjij are the sources of h, u\ and pij, 
respectively, and i,j = l,---,d. The explicit expression for 6 is given by, 

{d d \ 

-iX [/i(x, t) - h(t)] - i ^2 AW ~ i ^2 VijPij ( • (36) 
i=i i<j=i J 

We also define qijk by, qijk = h XiXjXk , which will be used later. Considering the zero-surface tension 
limit of the KPZ equation. The time evolution of the height /i(x, t) and its derivatives are given by, 

d 

d 

Ui,t = a^2uip u + f Xi , (38) 
l=i 

d d 

Pij,t = « XI PuPlj + « u Mv + fxw ■ (39) 
1=1 1=1 

Using Novikov's theorem 18 ' 19 (see Appendix D) we can write down the following identities, 

d 

</e> = -iXk(o)z - i mik"(o)z , (40) 
i=i 

(f Xi @) = -mk"{o)z , (4i) 

d 

if XiXi @) = -iXk"(0)Z - i mk""{0)Z , (42) 

1=1 

(fx i x j @} = -iVijk""(0)Z i^j, (43) 
where, for example, in d = 3 we have 

fc(x - x') = 2D D(x - x) , 

m = m = (^72 > 

fc'(o) = mo) = M°) = Mo) = o , 

k = k XX (0) = kyy(0) = k zz (0) 



-4D 
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If we write down the above equations for {d + 1) dimensions, the only change would be having d 
components for the arguments appearing in k, k", and k"" . 

Differentiating the generating function Z with respect to t, and using Eqs. (|37|) - (43) and the 
following identity, 

iZ Xl - i\Z m - i E V-i z mi = E VijiVijl®) > ( 44 ) 

i i<j 

we find that the time evolution of Z is governed by 

Z t = iXj(t)Z Z »ih ~ E Z m Viflnunu ~ ^K®) 2 + E $k"(0)Z 

l l l,i<j I 

-2\J2mik'\o)z- I E to + E^I k ""^) z > ( 45 ) 

I \l,k Kk j 

where, j(t) = hf. If we Fourier transform Z with respect to A, fa, and rjij, we obtain the governing 
equation for the joint PDF of h, U{, and pij, 

dX -p-r dm i — r drjij 



2vr xx 2vr 2vr 

x exp I iX [h(x, t) + h(t)] + i E + i E r^fte f Z(X, (ii,rnj,Xi, t) . (46) 

[ l l<k ) 

Using Eqs. (|4"5|) and (j4*6|) . the governing equation for the time evolution of P(h,Ui,Pij,t) in d spatial 
dimensions is then given by, 



l l Lk<m I 



+ 2*"(0) E + *""(0) E Ppunu ~ 2fc,W (°) E P PuP kk ■ (47) 

/ Z<fc Kk 

Equation ()47p enables us to obtain the governing equations for the time evolution of the moments of 
the height and its derivatives. The resulting equation has the following general form, 

^(h^AB) = -n 0l (t)(h^AB)-^Y,Ch n ^ 1 ABuf)+a E n km (h^AB P -^) 

01 1 l l,k<m Pkm 



h n °AB ^ ,h n °AB< 



-a 



E(^ ^/) + ^(0)n (no-l)(/i no - 2 A J B)-A : ,/ (0)E^(^-l)( — ) + 2k"(0) E n n u ( 



n 2 

I I "l I 



U l I Pll 



h n o A R h n ° A R 

+ k""(0) E n lk (n lk - 1)(^4^) + 2k"" (0) E n u n kk (^^-) , (48) 
i<k Plk i <k P"P kk 
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where 

i=i 

11 ll/v ; - 

i<j 

By using various values no, raj, and n^, coupled differential equations that govern the evolution of the 
moments are constructed. 
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Appendix B 



In this Appendix we prove the identity, (pij exp (—i\h — ij^i = 0. As shown in Appendix 

C, this identity is crucial to deriving Eq. (20). We also determine the height moments exactly and 
show that, {h XiXj Q} = for i ^ j, by considering an initially flat surface. However, it can be shown 
that this identity also holds for the general (d + 1)— dimensional surfaces. Setting no = ni = riij = 0, 
we obtain from Eq. (|48|) . 

-a£>H> = -«<V-u)=0=*(V-u)=0. (49) 
i 

The main aim is to calculate the moments {ptj exp(—i\h — ^i u i))- I n order to do so, we must 
follow several steps. First, we must determine such moments as, (pijUiUj) for i ^ j. Using Eq. (08]), 
we have 

d 

jr-iuiUjPij) = a^iuiUjPnpij) - ^^{mujpupij) = a^L^iU^pupij - puPij)) . (50) 
m I i i 

Inspecting the right-hand side of Eq. (|50p . we see that the terms with I = i ox I = j cancel one 
another. Noting that, for now, we have restricted our attention to (3 + 1)— dimensions, Eq. (|50|) is 
written as 

d 

TjjUi'ijP;;) = {UiUjipiiPij - PuPij)) i^j^l. (51) 
Inserting the right-hand side of Eq. (|5ip in Eq. (|48p . one finds 

■^{uiUjipupij - puPij)) = a^iuiUjPkiPkipij) + a^2{uiUjPkjPkiPu) 

k k 

- a^2(uiUjp kk piipij) - a^2(uiUjp k iPkiPij) - a ^{uiUjpkiPkjPu) + u^2(uiUjp kk PijPu) . (52) 

k k k k 

It can then be seen that for i ^ j ^ I the right-hand side of Eq. (|52p is zero which, when utilized, 
yields the following equation for a flat initial surface, 

(uiUjipupij - puPij)) = , (53) 

and 

(uiUjPij) = . (54) 

It can also be shown by induction that all the moments {h n °p^ u^vl^ 1 ) are identically zero. Therefore, 
we conclude that, (p^ exp (—i\h — iJ2i ^i u i)) = 0- 
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Appendix C 

Using the identities that we derived in Appendix B, we derive the joint PDF for the KPZ equation 
in the limit of zero surface tension. For clarity, we restrict ourselves to the (3 + 1)— dimensional 
surface, but generalization of the results to any (d + 1)— dimensional surface is quite straightforward. 
The zero-surface tension limit of the KPZ equation in (3 + 1)— dimensions has the following form (for 
more clarity, we replace the h(x,t) with the new notation h(x,y,z,t) in (3 + 1)— dimensions), 

h t {x,y,z,t)--a{h 2 x + h 2 y + h 2 z ) = f. (55) 

Letting, 

h x = u, h y = v, h z = w (56) 
and differentiating the KPZ equation with respect to x, y, and z, we have 

u t = a(uu x + vv x + ww x ) + f x , 

V t = 6i(uUy + VVy + WWy) + fy , 

Wt = a(uu z + vv z + ww z ) + f z , (57) 

for the height h and the corresponding velocity fields. The generating function Z(X, fj,±, ^2, x, y, z, t) 
is defined so as to generate the height and velocity field moments. By defining 6 as 

= exp {— i\[h(x, y, z, t) — h(t)] — ifi\u — i^v — i[i-$w} , (58) 

the generating function is written as, Z(X, ^2, A*3, x, y, z, t) = (0). Using the KPZ equation and its 
derivatives, the time evolution of Z(X, fj,2, fJ-3, x, y, z, t) is then written as, 

Zt = i^f(t)XZ — -iXa((u 2 + v 2 + w 2 )@) — iafii((uu x + vv x + ww x )Q) — ia^duiiy + vv y + ww y )Q) 

- iafi 3 {(uu z + vv z + ww z )Q) - iX(fO) - im{f x &) - i^ify®) - V 3 (/ 2 6) , (59) 
where j(t) = h t = ^a(u 2 + v 2 + w 2 }. Utilizing the statistical homogeneity of the quantities, we obtain 

Z x = {(-iXu - ifiiu x - ifi 2 v x ~ i^3Wx)Q) , (60) 
z y = {(—i^v - ifiiUy - i/J,2V y - ifi 3 w y )Q) , (61) 
Z z = ((—iXw - ifi\u z - ifi2V z - i[J,3W z )@) . (62) 

Because we analyze the system only in a time regime in which the cusp singularities (sharp valleys) 

have not yet formed, the order of the partial derivatives can be exchanged, 

d 2 h _ d 2 h 
dxidxj dxjdxi ' 
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and therefore, v x = u y , w x = u z , and w y = v z . Keeping the definition of 0, Eq. (|58|) . in mind, we may 
write 

d 

i- — {{-iniu x - ifJ, 2 v x - ifjL 3 w x )@} = (u x @) - im(uu x &) - i[j, 2 (uv x @) - i/j, 3 (uw x @) . (63) 



Using Eqs. (|60|) and (|63|) we obtain, 

d 

(u x Q) - im(uu x @) - ifi 2 (uv x Q} - i^{uw x Q) = -A- — (uQ) = -iXZ^^ (64) 
In a similar manner we find for the y— and z— directions that, 

- iXZ^^ = (v y Q) - ifj,i(vu y @) - i(x 2 (vv y @) - i/i 3 (vw x Q) , (65) 

- iXZ^^ = {w z @) - ifj,x(wu z @) - ifi 2 {wv z @) - i(j, 3 (ww z Q) . (66) 

Using Novikov's theorem 18,19 (see Appendix D) the expressions for (/0) and (f Xi Q) are written in 
terms of Z (see Appendix A). Therefore, we obtain 

Z t = ij(t)XZ - -iXaZ^^ - -iXaZ - ^iXaZ^^ 

- a(u x e) - a(v y Q) - a(w z Q) - X 2 k(0, Q)Z + MiMO, 0)Z + fi 2 2 k xx (0, 0)Z + f4k xx (0, 0)Z . (67) 
The {u Xi Q) terms are written as 

i i 

(u x Q) = — (@) x H ((iXu + ip, 2 v x + ifi 3 w x )Q) 

= -i—Zn - —(v x Q) - ^K6) , (68) 
and, similarly, for (v y Q) and (w z Q) we have 

( Vy Q) = -i—Z m - ^(u y @) - ^(w y Q) , (69) 

fJ-2 fJ-2 

{WzQ) = -i±z^ - ^(u z e) - ^( Vz e) . (70) 

The terms {h XiXj Q) (i ^ j), that appear in the equation for the time evolution of Z (for example, 
(u y Q)), prevent us from writing the Z— equation in a closed form. Fortunately, as was shown in 
Appendix B, such terms are zero for a flat initial condition. Therefore, the generating function Z 
satisfies the following equation 



(t)XZ - ^iXaZ^ lfll - ^iXaZ^ 2IM2 - ^iXaZ^^ + ia^-Z^ - ia-^-. 
A 



ia—Z„ 3 - X 2 k(0, 0)Z + n\k xx {0, 0)Z + /£k xx (0, 0)Z + f4k xx (0, 0)Z . (71) 
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We now solve Eq. (fTTj) assuming a flat initial condition, h(x, y, z, 0) = u(x, y, z, 0) = v(x, y, z, 0) 
w(x, y, z, 0) = 0, which is equivalent to writing 



P(h,u,v,0) = S(h)S(u)6(v)S(w) 



(72) 



which means that 

Z(0,0,*) = 1. (73) 
An efficient way of solving Eq. (|7ip is to factorize Z in the following manner, 11,15,16 

Z(A,Mi,M2,A»3,f) = J Fi(A,//i,t)F 2 (A,/i 2 ,t)i ? 3(A,^3,i)exp[-A 2 fc(0)t] . (74) 

Then, by inserting Eq. (|79p in Eq. (|7ip we obtain 

F U F 2 F 3 + F x F 2t F 3 + FiF 2 F 3t = i 1 {t)\F l F 2 F z - -iXaF 2 F lfJiltil F 2 F 3 - -iXaF^^Fs 



Ml 



Therefore, 



^2 ^3 

+ M"(0)FiF 2J F2 + /^"(O^i^ 



-iA<xF„„ + ia-i^ + [/i 2 /c"(0) - iaA/c"(0)t]F 
2 n 



(75) 



(76) 



with the initial condition, F(X,fi,0) = 1. This indicates that the height gradients in the three di- 
mensions evolve independently of one another before the emergence of cusp singularities, and that 
they are only coupled with the height field. By Fourier transforming Eq. (|76p with respect to ^ a 
simpler partial differential equation of order one is obtained, which can be solved by the method of 
characteristics. 11 The solution of Fj is then given by 



Fj(n, A, t) = < 1 — tanh 



tJ2ik xx (0)aX 



" 1/4 I 1 ,/// w 9 1 9 l2ik xx (0) 

exp { iak"(0)Xt 2 - -i[i 2 \\ -^r- 2 - tanh 



qA 



t^/2ik xx (0)aX 
(77) 



Therefore, we obtain the generating function in (3 + 1)— dimensions, 

-3/4 | | 



Z(X, fii, // 2 , t) = < 1 — tanh 



tJ2ik xx (Q)aX 



exp \ ™t(Mi + ^2 + ^sjy — ^ — 



tanh 



tJ2ik xx (0)aX 



x exp 



-^aA;"(0)Ai 2 - k(0)X 2 t 



(78) 



By inverse Fourier transformation of the generating function Z, the PDF of the height fluctuations, 
P(h,u,v,w,t), is then determined. 
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As mentioned earlier, generalizing the above results to any (d + 1)— dimensional surface is straight- 
forward. For example, for d = 2 all we need to do is setting /X3 = in Eq. (178 j) . For a general d we 
have the following expression for Z: 

Z(X, m, ■ ■■ ,Hd,t) = Fi(\,fj,i,t) • • • F d (X,fi d ,t) exp[-A 2 A;(0)t] , (79) 

where Fa is defined as before. 
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Appendix D 

In this Appendix we prove a simple version of the theorem by Novikov 18 ' 19 that we use in this 
paper. As before, we define the generating, Z(X, fa, X\, t) = (0(A, fa, X{, t)), for the fields h = h(x, t) — h 
and Ui = h Xi = dh/dxi by 



6 = exp | -i\[h(x, y, z, t) - h(t)] - i ^ faUi j , 



(80) 



where, as usual, A and fa are the sources of h and fa, respectively, and i, j = 1, • • • ,d. Starting with 
the equation, 

d t h{x,t) = -C[h(x,t)] + f(x,t) (81) 

where /(x, t) is a term representing the noise with Gaussian correlations in space and white noise in 
time: 

(/(x, t)/(x', t')) = 2D D(x - x')5(t - t') , (82) 

we can write, 

<F[/]/(x,t)> = J dx'dt' (/(x,i)/(x' 5 t')> ( gTpfk ) • ( 83 ) 



If we take F = O, we obtain 

d 



dF [f] _ ■ \ dh(x,t) pm x ^ dfa(x,t) pm 



(84) 



Using Eq. (81) we can write, 

r *£[/i(x,t)]dt'+ /' 
which, when differentiated with respect to /, gives, 



h(x, t) = h(x, t ) - f £[h(x, t))dt' + /' /(x, t')dt' (85) 

Jtn Jtn 



S/(x,t) Jto /(x,t) 

Using Eqs. (84) and (86) in the limit, t — t' — > 0, invoking casuality, and considering the same 
expressions for, J2i=i ^idui(x,t)/df(x' ,t')F[f], we obtain 

d 

(F[f]f{x,t)) = -2i\D D(0){F[f])-2iJ2^D D"(x-x')(F[f]) . (87) 

i=l 

Using the last expressions and k(x) = 2DqD(~x), we obtain Eq. (|4(jp . 
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Figure 1: Different time snapshots of the velocity fields of randomly driven one-dimensional burgers 
equation 



Figure 2: Plot of f + vs a for the KPZ equation in the strong coupling and before the emergence of 
sharp valleys for various times, t/t* in (2 + 1)— dimensions (color online). 
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Figure 3: Logarithmic plot of iV+ vs t for the KPZ equation in the strong coupling, and before the 
emergence of sharp valleys in (2 + 1)— dimensions (color online). 
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